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Abstract 

Zr and Sc precipitate in aluminum alloys to form the Al3Zr2.Sci_a; compound which, for low 
supersaturations of the solid solution, exhibits the LI2 structure. The aim of the present 
study is to model at an atomic scale the kinetics of precipitation and to build mesoscopic 
models so as to extend the range of supersaturations and annealing times that can be sim- 
ulated up to values of practical interest. In this purpose, we use some ab initio calculations 
and experimental data to fit an Ising type model describing thermodynamics of the Al-Zr-Sc 
system. Kinetics of precipitation are studied with a kinetic Monte Carlo algorithm based on 
an atom-vacancy exchange mechanism. Cluster dynamics is then used to model at a meso- 
scopic scale all the different stages of homogeneous precipitation in the two binary Al-Zr and 
Al-Sc alloys. This technique correctly manages to reproduce both the kinetics of precipi- 
tation simulated with kinetic Monte Carlo as well as experimental observations. Focusing 
on the nucleation stage, it is shown that classical theory well applies as long as the short 
range order tendency of the system is considered. This allows us to propose an extension of 
classical nucleation theory for the ternary Al-Zr-Sc alloy. 

Introduction 

Transition elements like Mn, Cr, Zr or Sc are usually added to aluminum alloys so as to form 
small ordered precipitates which increase the tensile strength and inhibit recrystallization. 
As these properties directly depend on the density and the size of the particles, an accurate 
control of the precipitation process is necessary so as to optimize the material. In particular, 
one can use a combined addition of solute elements so as to change the precipitate distribution 
and thus improve the alloy properties. Such an example for aluminum alloy is the addition 
of Zr and Sc elements. Both elements, when added separately, increase the tensile strength 
and the recrystallization resistance, but the effect is even better for the combined addition 
as there are more precipitates which are smaller and less sensitive to coarsening [1-6]. So as 
to better understand the precipitation kinetics in Al-Zr-Sc we develop a multiscale approach 
combining atomic and mesoscopic simulations. 

At the atomic scale, kinetic Monte Carlo (KMC) simulations are the suitable tool to 
study precipitation kinetics. Thanks to a precise description of the physical phenomenon 
leading to the evolution of the alloy, such simulations allow to predict the kinetic pathways 
in full details. These predictions are valuable for a ternary alloy like Al-Zr-Sc because few 
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information on the precipitation kinetics is available whereas the diversity of the possible 
kinetic pathways is richer than for a binary alloy. Nevertheless, one drawback of this atomic 
approach is the needed computational time: only the early stages of precipitation and the 
high supersaturations can be simulated. There is thus a need to develop mesoscopic models 
so as to study alloys and heat treatments compatible with technological requirements and 
to allow a comparison between simulated kinetics and experimental data. These mesoscopic 
models, which include cluster dynamics (CD) and the ones based on classical nucleation 
theory (CNT), describe the alloy by the means of the cluster size distribution. Their input 
parameters are not always directly available but can be obtained from the atomic model. 
Their main drawback is that they are not as predictive as atomic models, especially for 
ternary alloys where one has to know the kinetic pathway before modeling the precipitation 
kinetics. Therefore, a multiscale approach combining atomic and mesoscopic simulations 
sounds valuable. It allows to use KMC simulations to obtain the missing information required 
by mesoscopic models and thus extend the range of supersaturations and heat treatments 
that can be studied. 

Kinetic Monte Carlo (KMC) 

Atomic model 

Al-Zr-Sc thermodynamics. The system we want to model is the ternary Al-Zr-Sc alloy. The 
aluminum solid solution has a face-centered-cubic (fee) structure and the forming precipitates 
are AlaZr, AI3SC, and AI^ZyxSci-x, all precipitates having the LI2 structure.^ This structure 
relies on an fee underlying lattice. Therefore a rigid lattice model can be used so as to 
describe the different configurations of the system. This model is well suited to the Al-Zr-Sc 
alloy because the lattice parameters of the forming precipitates are really close to that of 
the aluminum matrix and precipitates remain coherent at least for diameters smaller than 
10 nm. Therefore elastic relaxation as well as loss of coherency can be neglected in a study 
focusing on the first stages of the precipitation kinetics. 

Atoms are thus constrained to lie on an fee lattice configurations of which are described 
by the occupation numbers p^. = 1 if the site n is occupied by an atom of type i and 

= if not. Energies of such configurations are given by an Ising model with first and 
second nearest neighbor interactions. Thus, in our model, the energy per site of a given 
configuration is 

^''^s n,m ^-'^s r,s 

where the first and second sums respectively runs on all first and second nearest neighbor 
pairs of sites, Ng is the number of lattice sites, e\f and e^f are the respective effective energies 
of a first and second nearest neighbor pair in the configuration 

Such a rigid lattice model has been previously developed for Al-Zr and Al-Sc binary 
alloys [12]. First and second nearest neighbor interactions have been fitted so as to reproduce 
the free energies of formation of AlsZr and AI3SC compounds in the LI2 structure as well as 
Zr and Sc solubility limits in aluminum.^ So as to study precipitation in the ternary Al-Zr-Sc 

^AlaZr and AlsZr^^Sci-a; for 0.8 < x < 1 have the stable DO23 structure [1,7,8]. Nevertheless for small 
supersaturations of the solid solution, precipitates with the LI2 structure nucleate and grow first. Precipitates 
with the DO23 structure only appear for prolonged heat treatment and high enough supersaturations [9-11]. 

^For Al-Zr interaction, as we want to model precipitation of the metastable LI2 structure, we use the 
metastable solubility limit that we previously obtained from ab initio calculations [13]. 



2 



system, this atomic model needs to be generalized by including interactions between Zr and 
Sc atoms. As no experimental information is available, we use ab initio calculations to deduce 
the corresponding parameters. The formation energies of 19 ordered compounds in the Al- 
Zr-Sc system are calculated with the full potential linear muffin tin orbital method [14] in the 
generalized gradient approximation [15]. We then use the inverse Connoly- Williams method 
[16] to obtain the two missing parameters of the atomic model, i.e. e^rSc ^ZrSc fitted on 
this database of formation energies. We obtain an order energy cj^rSc — ^zrSc"~|^ZrZr~|^ScSc ~ 
237 meV corresponding to a strong repulsion when Zr and Sc atoms are first nearest neighbors 
and a slight attraction when they are second nearest neighbors (c^z/sc ~ —2.77 meV). With 
such interactions, an ordered ternary compound AlgZrSc is stable at OK. It is based on 
the LI2 structure with Zr and Sc atoms being ordered on the minority sublattice so as 
to form only attractive second nearest neighbor ZrSc interactions. When increasing the 
temperature, this structure partially disorders and leads to an AlaZr^-Sci.^ compound where 
< X < 1 is a variable quantity. This compound has a LI2 structure too and the atoms on 
the minority sublattice can equally be Zr or Sc. The order-disorder transition temperature 
estimated within Bragg- Williams approximation is T ~ 70K. Thus, at the temperatures we 
are interested in, i.e. above the ambient temperature, our atomic model predicts that only 
the AlaZr^-Sci-^. compound is stable. This agrees with experimental observations [1, 8, 17] 
performed with transmission electron microscopy (TEM) which show that precipitates in the 
ternary system have the structure described above. 



Al-Zr-Sc kinetics. We introduce in the Ising model atom-vacancy interactions for first near- 
est neighbors: e^iV' ^ZrV' ^'^'^ ^ScV deduced from vacancy formation energy and binding 
energies with solute atoms in aluminum [12]. 

Diffusion can then be described through vacancy jumps. The vacancy exchange frequency 
with one of its twelve first nearest neighbors of type a is given by 

— w)' ^'^^ 

where Va is an attempt frequency and the activation energy E*^^ is the energy change required 
to move the a atom from its initial stable position to the saddle point position. It is computed 
as the difference between the contribution e^^ of the jumping atom to the saddle point energy 
and the contribution of the vacancy and of the jumping atom to the initial energy of the 
stable position. This last term is obtained by considering all bonds which are broken during 
the jump. The six kinetic parameters and e^^ are fitted so as to reproduce Al self diffusion 
coefficient and Zr and Sc impurity diffusion coefficients [12]. 



Atomic simulations 

Kinetic Monte Carlo algorithm. We use time residence algorithm to run KMC simulations. 
The simulation boxes contain Ns = 100^ or 200^ lattice sites and a vacancy occupies one of 
these sites. At each step, the vacancy can exchange with one of its twelve first nearest neigh- 
bors, the probability of each jump being given by Eq. (j2)). The time increment corresponding 
to this event is 

_ 1 1 

'-iV.(l-13xO)Cv(Al)^;^l^r._^' 

where = x^j- + is the solute nominal concentration of the simulation box and Cv(Al) 
the real vacancy concentration in pure Al as deduced from energy parameters. 
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Criteria used to discriminate atoms belonging to the solid solution from those in LI2 
precipitates are the same as the ones developed in our study of the precipitation kinetics 
in the two binary Al-Zr and Al-Sc alloys: perfect order is assumed for LI2 cluster'^ and a 
critical size is defined so as to differentiate unstable clusters from precipitates. 

Precipitation kinetics in the ternary alloy. Using this atomic diffusion model, we run KMC 
simulations of the annealing of supersaturated Al-Zr-Sc sohd solutions. Precipitates appear- 
ing in the simulation box are AlsZr^-Sci-a; ternary compounds having the partially ordered 
LI2 structure predicted by thermodynamics. These precipitates are inhomogeneous (Fig. 
their core is richer in Sc than in Zr whereas this is the opposite for the external shells. 
Zr concentration slightly decreases away from the core and the intermediate shells have 
almost the AI3SC stoichiometry. As for the external shells, they are strongly enriched in 
Zr. This structure predicted by our atomic model agrees with experimental observations 
recently made with high resolution electron microscopy [18, 19] and with three-dimensional 
atom probe [20,21]. 




r / a 

Figure 1: Precipitate obtained after the annealing during 1.2 s at 550°C of an aluminum 
solid solution containing 0.1 at.% Zr and 0.5 at.% Sc simulated with KMC (for clarity, only 
Zr and Sc atoms are shown). The corresponding radial concentration profile shows the Zr 
enrichment of the periphery compared to the core of the precipitate. 

Following the composition and the size of all clusters present in the simulation box 
(Fig. 12)), we see that small clusters, which should be sub-critical, can have all compositions 
between the two AlsZr and AI3SC stoichiometric compounds. As for over-critical clusters, 
the smaller ones are richer in Zr than the larger ones. These small clusters, which have 
just nucleated, grow by absorbing Sc atoms, leading to the tail of the distribution in the Sc 
richer part. As there is no diffusion inside precipitates (because of a high vacancy energy 
of formation and of the energetic cost associated with the creation of antisite defects), Zr 
and Sc atoms cannot homogenize. Therefore this distribution is consistent with the previous 
observation that the core of the precipitates can contain Zr whereas the intermediate shells 
contain almost exclusively Sc. This shows that Zr plays a role during nucleation but that 
growth is controlled by Sc, because this element diffuses much faster than Zr in the matrix. 
As precipitation is going on, the solid solution depletes in Sc. Therefore clusters can only 
grow by absorbing Zr atoms, leading to a move of the whole distribution towards the Zr 

^Zr and Sc atoms belong to a LI2 cluster if all their first nearest neighbors are Al atoms and at least one 
of their six nearest neighbors is a Zr or Sc atom. 



4 




Figure 2: Projection in the plane xzr, nx of the precipitate distribution in a sohd solution 
of composition x^^ = 0.5 at.% and Xg^ = 0.5 at.% annealed at T = 550°C. 

richer part (Fig. This is consistent with the fact that the external shells of the precipi- 
tates are richer in Zr. The precipitate growth by Zr absorption occurring in a second stage 
leads to a slowdown of the precipitation kinetics. 

In agreement with TEM observations [1,4], these atomic simulations show that a Zr 
addition to an Al-Sc alloy leads to a higher density of precipitates, these precipitates being 
smaller (FiglH)). Our atomic model allows to understand this effect. A Zr addition increases 
the nucleation driving force and thus the number of precipitates. As these ones mainly grow 
by absorbing Sc, one obtains smaller precipitates because the number of growing clusters 
is higher and there is less Sc available for growth at the end of the nucleation stage. The 
precipitate growth by Zr absorption is really slow but plays an important role too as it 
leads to the formation of an external Zr enriched shell which is responsible for the good 
resistance against coarsening. This kinetic effect due to the addition of an impurity having 
an attractive interaction with solute atoms (c^zj-sc ^ 0) been already noticed by Soisson 
and Martin [22]. 

It is interesting to notice that a Sc addition to an Al-Zr alloy has quite distinct effects 
(Fig. 12)). As Sc atoms plays a key role for nucleation and growth, such an addition leads to 
an increase of the precipitate density as well as to an increase of the precipitate size. Indeed 
the nucleation stage is shorter and the number of nucleating precipitates is higher. But this 
increase of the nuclei density does not consume the whole Sc in the solid solution and at the 



5 




Figure 3: Kinetics of precipitation at T = 450°C: variation with the Zr and Sc concentration 
of the number Np of precipitates and of their mean size < nx >■ The critical size to 
discriminate precipitates from sub-critical clusters is — 6. 

end of the nuclcation stage, enough Sc remains so as to accelerate precipitate growth and 
coarsening, thus leading to bigger precipitates. In this case, one has to consider not only the 
thermodynamic contribution of the impurity addition, but its kinetic contribution too so as 
to conclude on its effects on the precipitation kinetics. 

Cluster dynamics (CD) 

KMC simulations allow us to understand precipitation kinetics in the ternary Al-Zr-Sc alloy 
in great details. Nevertheless, these simulations are restricted to short annealing times and 
high enough supersaturations, preventing any comparison with experimental data. Starting 
from the same atomic model, we build a mesoscopic modeling [23] which is based on a set of 
rate equations describing in the same framework the three different stages of precipitation, 
i.e. nucleation, growth and coarsening. It requires only a limited number of parameters, the 
interface free energy and the solute diffusion coefficients, which can be quite easily deduced 
from atomic parameters. The range of supersaturations and annealing times that can be 
simulated are thus extended for both binary Al-Zr and Al-Sc systems allowing a comparison 
with experimental data. 

Mesoscopic model 

In its strict sense, CD rests on the description of the alloy undergoing phase separation 
as a gas of solute clusters which exchange solute atoms by single atom diffusion [24-26]. 
Clusters are assumed to be spherical and are described by a single parameter, their size or 
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the number nx of solute atoms X they contain. In such a description, there is no precisely 
defined distinction between the solid solution on the one hand and the precipitates on the 
other hand, at variance with the CNT: the distribution of cluster sizes is the only quantity 
of interest (for a detailed discussion see Ref. 27). 



Master equation. The CD technique describes the precipitation kinetics thanks to a mas- 
ter equation giving the time evolution of the cluster size distribution [24-26]. When only 
monomers can migrate, which is the case for Al-Zr as well as Al-Sc systems [12], the prob- 
ability to observe a cluster containing nx solute atoms obeys the differential equations 



dC 

~ Jnx-l^nx ~ Jnx^nx+ly V > 2 (4a) 

dCi 

= -2Jl^2 - ^ J„x-nx+l, (4b) 
nx>2 

where the fiux J^^^nx+i from the class of size nx to the class nx -l- 1 is written 

^nx~*«x+l Pnx^nx '-^^nx+l^^x+l ' (^) 

Pnx being the probabihty per unit time for one solute atom to impinge on a cluster of size 
nx and for one atom to leave a cluster of size nx- 



Condensation rate. When the solute long-range diffusion controls the precipitation kinetics, 
the condensation rate is obtained by solving the diffusion problem in the solid solution 
around a spherical precipitate. For a cluster of radius r^^, this condensation rate takes the 
form [28, 29] 

Dx 

Pnx = ^^r^^—Ci, (6) 

where Dx is the diffusion coefficient of the solute at infinite dilution of the solid solution and 
Q is the mean atomic volume corresponding to one lattice site. 



Evaporation rate. The evaporation rate is obtained assuming that it is an intrinsic property 
of the cluster and therefore does not depend on the solid solution surrounding the cluster. 
This means that the cluster has enough time to explore all its configurations between the 
arrival and the departure of a solute atom. Thus should not depend on the nominal 
concentration of the solid solution and can be obtained by considering any undersaturated 
solid solution of nominal concentration x^. Such a solid solution is stable. Then there should 
be no energy dissipation. This involves that all fiuxes Jnx^nx+i equal zero. Using equation 
El one obtains 

r< (^0 \ 

/0\ Jd / 0\ ^^''^x \ X/ f ''7\ 

L'nx+ll^^xJ 

where overlined quantities are evaluated in the solid solution at equihbrium. This finally 
leads to 

Dx 

anx+i = 47rr„x-^exp [(G„x+i - -Gi)/kT], (8) 

where Gn^ is the free energy of a cluster containing nx atoms. It can be divided into a volume 
and an interface contributions, which involves the following expression of the evaporation 
rate: 

anx+i = 47rr„x^ exp [{SGnY^^a^ {{nx + l)'/V„x+i - nx'^V^x - ^i) /kT] , (9) 
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where (T„x is the interface free energy of a cluster containing nx solute atoms and a the 
lattice parameter.^ 

Looking at the expression IHl of the condensation rate and the expression El of the evapo- 
ration rate, one sees that the only parameters needed by CD are the diffusion coefficient and 
the interface free energy. There is no need to know the nucleation free energy in opposition 
to other mesoscopic models based on CNT. In CD, thermodynamics of the solid solution is 
described thanks to a lattice gas model and therefore the nucleation free energy results from 
this description and is not an input of the modeling. 

Interface free energy. The interface free energy can be deduced from the atomic model [12] 
by computing within the Bragg- Williams approximation the free energies corresponding to 
planar interfaces between the aluminum solid solution and the LI2 precipitates for the three 
most dense packing orientations {111}, {110} and {100}. We obtain ctioo < cthq < o'lu, 
indicating that precipitates mainly show facets in the {100} direction and that facets in 
the {110} and {111} directions are small. As the difference between aioo, o"iio and am is 
decreasing with temperature, precipitates are becoming more isotropic at higher tempera- 
tures. This is in agreement with the shapes of the precipitates observed during the atomic 
simulations, as well as with experimental observations done by Marquis and Seidman [30]. 

So as to get the input parameters needed by CD, one has to deduce an isotropic average 
interface free energy a from these planar interface free energies. This can be done using 
the Wulff construction [31,32]: a is defined so as to give the same interface free energy for 
a spherical precipitate having the same volume as the real faceted one.^ For temperatures 
ranging between K and the aluminum melting temperature (933.5 K), a is varying be- 
tween 124 and 92 mJ.m~^ for AlsZr precipitates and between 138 and 95 mJ.m~^ for AI3SC 
precipitates. 

The interface free energy a does not depend on the size of the cluster and corresponds to 
the asymptotic limit of the quantity used to define the evaporation rate (Eq. IH)). A direct 
calculation considering the partition functions of small clusters [12] shows that slightly 
deviates from this asymptotic value at small sizes (Fig.|3)). Nevertheless, CD is very sensitive 
to this input parameter because its thermodynamic description only relies on it. Therefore, 
the size dependence of cr„x has to be considered. In this purpose, for clusters containing 
no more than 9 solute atoms, the interface free energy is computed precisely using Eq. 20 
of Ref. 12 and for clusters of size nx > 10, the interface free energy is obtained using an 
extension of the capillary approximation. 



where c and d respectively correspond to the line and point contributions [33]. The asymp- 
totic value a used is the one previously calculated whereas coefficients c and d are obtained 
by a least square fit of the exact expression for sizes 5 > nx > 9. 

Precipitation kinetics 

In order to compare precipitation kinetics obtained from KMC and from CD simulations, an 
arbitrary threshold size is used to define precipitates so as to be able to calculate mean 
quantities like the precipitate size and density. Below this threshold size, clusters are assumed 

more detailed description of the derivation of Eq. ^is given in Ref. 23. 
^Details of the calculation can be found in the appendix A of Ref. 12. 




(10) 
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Figure 4: Variation with the cluster size nsc of the interface free energy between the sohd 
solution and AI3SC at three different temperatures (300, 400, and 500°C). Symbols correspond 
to cr„g^ as given by the direct calculation of the cluster formation free energy (Eq. 20 of Ref. 
12) and lines to the capillary approximation (Eq. ITUI) as well as to its asymptotic limit. 
Dotted lines indicate the range of sizes where the coefficients c and d have been fitted. 

to be invisible. When comparing with experimental data, we choose = 27 corresponding 
to a precipitate diameter equal to 1.5 nm which is supposed to be the smallest size which 
can be discriminated by TEM. This threshold size can thus be different from the critical size 
of CNT: the only important thing is to use the same criterion to define precipitates when 
comparing atomic and mesoscopic simulations. 

Comparing the precipitation kinetics obtained with KMC and CD simulations, one sees 
that CD manages to reproduce the variations of the precipitate density and of their mean 
size (Fig. EJ. A more thorough comparison [23] shows that CD catches variations with 
the supersaturations and the annealing temperatures of the precipitation kinetics. For low 
supersaturations, the agreement between both modeling is really good, whereas for high 
supersaturations, kinetics modeled by CD appear to be slightly too slow compared to the 
ones simulated with KMC. This delay only corresponds to a constant factor on the time scale 
and in all cases the prediction of the precipitate maximal density at the transition between 
the growth and the coarsening stages is correct. 

Precipitation kinetics simulated with CD can be compared too with the experimental data 
of Novotny and Ardell [34] and those of Marquis et al. [30,35,36] who studied an aluminum 
alloy having the composition x^^ = 0.18 at.%. For the three different temperatures T = 300, 
350 and 400°C, it appears that the CD equations manage to reproduce the variation with 
time of the mean precipitate radius (Fig. E)). Despite the fact that the time scales in real 
experiments are several order of magnitude larger than those in KMC, CD, with a single set 
of parameters, reproduces atomic simulations at short times and gives a safe extrapolation 
thereof to the range of annealing times that can be compared with experimental data. In 
Ref. 23, it was shown too that the precipitate size distributions simulated with CD agrees 
really well with the experimentally ones measured by Novotny and Ardell [34]. 

CD can be used so as to simulate the evolution of the electric resistivity during annealing 
kinetics. Some calculations of the electrical resistivity associated with a cluster distribution 
exist in the literature [39,40], but as a first approximation, one can consider that the resis- 
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Figure 5: Mean precipitate radius and precipitate density as a function of the aging time 
for a solid solution of composition Xg^ = 0.18 at.% at different temperatures (T = 300, 350, 
and 400°C) as given by CD and compared to experimental data [30, 34-36] and to KMC 
results [12]. The cutoff radius used for CD and KMC is r* ~ 0.75 nm {n^ = 27). 

tivity contribution p„ of a cluster of size n is proportional to its section, p„ = piv?^^. If one 
assumes that all clusters are contributing to the resistivity, the resistivity measured at time 
t in the phase separating system should be 

oo 

p(^) = pO^ + 5psc5^c„(^y/^ (11) 

n=l 

where the size distribution Cnif) is given by CD (Eq. HJ. is the resistivity of pure 
aluminum and is temperature dependent.^ The increase of resistivity with Sc content has 

^Measurements corresponding to experimental studies [37, 38] have been performed at different tempera- 
tures. When comparing with R0yset and Ryum experimental data [38], we use p^j = 27 nfim, whereas we 
use p^j = 28.4 nfim for the comparison with Zakharov data [37]. 
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Figure 6: Time evolution of the resistivity p experimentally observed [37, 38] and deduced 
from CD simulations for two solid solution of composition Xg^ = 0.12 and 0.24 at.%. 

been measured at 77 K by Fujikawa et al. [41,42] who give 5psc = 3400 nfim. Assuming 
that Matthienssen rule [43] is obeyed, this quantity does not depend on the temperature. 
When comparing resistivity predicted by CD with the experimental ones measured during 
precipitation kinetics [37,38] we obtain a good agreement (Fig. Ej), especially for the lowest 
temperatures: CD manages to reproduce the fast decrease of resistivity during the nucleation 
and growth stages as well as the slower variation which follows during the coarsening stage. 
For the highest temperatures (330°C for the solid solution of concentration Xg^ = 0.12 at.% 
and 400°C for Xg^, = 0.24 at.%), CD appears to be too fast compared to experimental data. 
This may arise from the approximation we are using for the contributions to resistivity of 
the various clusters as this size dependence may be too crude for small clusters. One should 
notice too that the evolution predicted by CD at these temperatures is really fast and that 
the resistivity drop appears at times (t ~ 10 s) which are too small to be precisely observed 
experimentally. 

CD is thus a powerful modeling technique to study precipitation in a binary alloy, allow- 
ing to link different scales. Starting from an atomic diffusion model, it gives quantitative 
predictions of precipitation kinetics in agreement with experimental data. It is possible to 
generalize this mesoscopic technique so as to study precipitation in a ternary alloy [44] but 
some work still needs to be done. In the case of the Al-Zr-Sc alloy, one of the main difficul- 
ties is to find a way for CD to describe the precipitate inhomogeneities. Therefore, in the 
following, we will not use this technique to model the whole kinetics of precipitation in the 
ternary alloy, but we will focus on the nucleation stage and see how CNT can be extended 
so as to predict the nucleation rate as well as the compositions x of the AlaZr^jSci-a; nuclei 
in this ternary alloy. 

Classical nucleation theory (CNT) 

CNT^ allows to predict the nucleation rate in a supersaturated binary solid solution. Com- 
parisons with KMC simulations [12,22,47-51] show that its predictions are correct as long as 
its input parameters are carefully evaluated. In particular, for the Al-Zr and Al-Sc alloys, it 
has been shown that one has to take into account the ordering tendency of the binary alloy 

^Detailed descriptions of the theory can be found in Ref. 29,45,46. 
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when calculating the nucleation driving force [12]. We will see how this can be done in a 
convenient way as well as how CNT applies to the two binary Al-Zr and Al-Sc alloys, before 
extending the theory to the ternary Al-Zr-Sc alloy. 



Binary Al-Zr and Al-Sc alloys 

Whereas CD describes the time evolution of the whole cluster size distribution (Eq. CNT 
assumes that the distribution is stationary in the solid solution, i.e. below a critical size n^, 
and does not provide any information on the distribution of larger clusters. The probability 
to observe a cluster containing nx solute atoms for nx < is thus given by 



= exp {-AGny,/kT), 



(12) 



where AG„^ is the formation free energy of the cluster.^ It can be evaluated thanks to the 
capillary approximation, dividing this formation free energy into a volume and a surface 
contribution, 



AG„,(4) 



Anx^G 



nuc/ 
Xv 



(367r 



,1/3 2/3 2 



(13) 



One thus sees that CNT needs one more input parameter than CD, the nucleation free energy 



[X 



XJ 



With this additional parameter and the assumption on the cluster size distribution in 
the metastable solid solution, CNT leads to the following prediction of the steady-state 
nucleation rate. 



J 



St 



1 d^AGn^ 




2nkT dnx^ 





exp 



AG 



rax* 



kT 



(14) 



where A^^ is the number of lattice site. The condensation rate takes the same expression^ 
as in CD (Eq. ^ and, like the formation free energy and its second derivative, it is evaluated 
for the critical size nX^ corresponding to the maximum of the formation free energy. 

Different thermodynamic approximations can be used to obtain the nucleation free en- 
ergy from the atomic model. Usually, one uses simple mean-field approximation, like the 
ideal or regular solid solution model which both rely on the Bragg- Williams approximation. 
Nevertheless, we showed in Ref. 12 that they lead to bad predictions of the CNT: the clus- 
ter size distribution in the solid solution corresponding to Eq. completely disagrees with 
that observed during KMC simulations. As the steady-state nucleation rate relies on this 
distribution, we observe the same discrepancy for the prediction of J'^* (Fig. [7j). So as to 
obtain a good agreement, one has to take into account order effects corresponding to the 
strong attraction existing between Al and Zr as well as Al and Sc atoms when they are first 
nearest neighbors and the strong repulsion when they are second nearest neighbors This can 
be done using the cluster variation method (CVM) [52, 53] to calculate the nucleation free 
energy. The cluster size distribution given by Eq. Elthen corresponds to the one observed in 
atomic simulations and CNT theory well reproduces the nucleation rate (Ref. 12 and Fig.[7j). 
Nevertheless, one drawback of CVM is that it does not lead to an analytical expression of 
the nucleation free energy. It is then not really convenient to combine it with a theory whose 



®The formation free energy can be related to the free energy Gn^ of the cluster by AG„x(a^x) = 
G„x — 2rixM(2;x)i where nix^) = {fJ-x{xx) ~ f^Aii^x)) /2 is the effective potential, i.e. a Lagrange mul- 
tiplier imposing that the nominal concentration of the solid solution is x^- 

^The monomer concentration Ci appearing in Eq. El has to be replaced by the nominal concentration 
in CNT. 
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Figure 7: Variation with the nominal concentration and the temperature of the steady-state 
nucleation rate J*** for AlsZr precipitation. Symbols correspond to KMC simulations and 
lines to CNT. On the left, different thermodynamic approximations are used to evaluate 
/\Qnuc. ^YiQ ideal and regular solid solution models, the cluster variation method (CVM) and 
the low temperature expansion (LTE) to second order (Eq. [T3j) . On the right, only LTE is 
used. 



main contribution is to give an analytical expression of the nucleation rate. Moreover, it 
would be interesting to understand thoroughly what is wrong with the ideal and the regular 
solid solution models and in particular why the regular model leads to worse predictions than 
the ideal one. In this purpose, we use low temperature expansions [54]. This thermodynamic 
approximation leads to analytical expression of the nucleation free energy and takes correctly 
into account the ordering tendency of the alloy. 

A low temperature expansion consists in developing the partition function of the system 
around a reference state, keeping in the series only the excited states of lowest energies. For 
an Al-Sc or an Al-Zr solid solution, this reference state corresponds to pure Al and the two 
first excited states are the transformation of an Al atom into an X solute atom and of a 
pair of Al second nearest neighbors to a pair of X atoms. When considering only the second 
order of the expansion, one can obtain analytical expressions of the different thermodynamic 
quantities in the canonic ensemble, thus depending on the nominal concentration of the 
solid solution. In particular, the nucleation free energy takes the form 



AG""'=(4) =kT [q(x°) -q{x'^)] + 3kT exp (2uj1}^/kT) 



q (x^)^ - q 



hT{\n [q(x°)]-ln[q(x^«)]} 



(15) 



where we have defined the function 



2r 

q (x) = (16) 



l + 24a;exp (2j^j^/kT) 



Using this expression of the nucleation free energy with CNT one obtains predictions of the 
nucleation rate as good as with CVM (Fig. [7j) compared to the rate observed during KMC 
simulations. In particular, the variations with the concentration of the solid solution and 
with the annealing temperatures are well reproduced. 

detailed description of the calculations can be found in Ref. 55. 
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So as to understand why simple mean-field approximations lead to bad approximations of 
the nucleation free energy, we develop the expression to first order in the concentrations 

^Gl-,{xl) ~ ^A:Tln (i^) + hrin + 6e-i^.'xATj (^o _ _ (^7) 

Doing the same development for the nucleation free energy calculated within the Bragg- 
Williams approximation (regular solid solution), we obtain 

AG-(4) - In (i^) + \kT In + (gJ^^ + Sc^x) (x^ - x^) . (18) 

Comparing Eq. El with Eq. we see that these two thermodynamic approximations de- 
viate from the ideal solid solution model by a distinct linear term. In the low temperature 
expansion (Eq. ITTj) . the nucleation free energy is only depending on the second nearest 
neighbor interaction and the coefficient in front of the concentration is positive. On the 
other hand, the Bragg- Williams approximation (Eq. ITS|) incorporates both first and second 
nearest neighbor interactions into a global parameter ujaix = ^^mx + '^^mx- This leads to 
a linear correction with a negative coefficient as ujaix < for both binary alloys. We thus 
demonstrate that the regular solid solution model leads to a wrong correction of the ideal 
model because it does not consider properly short range order. In the case of a LI2 ordered 
compound precipitating from a solid solution lying on a fee lattice, one cannot use such an 
approximation to calculate the nucleation free energy. On the other hand, Eq. is a good 
approximation and can be used to calculate the nucleation free energy even when the second 
nearest neighbor interaction cu^ix '^^^ known. Indeed, this parameter can be deduced from 
the solubility limit by inverting Eq. 4 of Ref. [12], leading to the relation 

41 = -hr In ix'^) + kT (^x^^' + fxt^''^ . (19) 

This relation combined with Eq. provides a powerful way for calculating the nucleation 
free energy from the solid solubility. 



Ternary Al-Zr-Sc alloy 

We now explore how CNT can be extended to the Al-Zr-Sc alloy. Going from the binary to the 
ternary alloy, the main difficulty is that the precipitates have the composition AlsZr^^Sci-^, 
where x is not known a priori. Therefore, the theory has to predict their composition. As 
we showed for the two binary Al-Zr and Al-Sc alloys that the capillary approximation is well 
suited to describe thermodynamics of the clusters, even the smallest ones. We use the same 
approximation for the ternary system. The formation free energy of a cluster AlsZr^Sci.^; 
containing nx solute atoms can then be written 

AG„,(xO„xO„x) = 4nxAG""^(x°„x°„x) + (Senf^nx'/^a'aix). (20) 

The nucleation free energy is depending not only on the composition of the supersaturated 
solid solution but also on the composition of the nucleus: 

AG-(x°„ xl, x) =\ {x [2^^l,{x) + /.^(x)] + (1 - x) [3/.^,(x) + ^^l,{x)] 

— X [3/i^'\(x2r, X^c) + l^zA^Zr^ 4c)] ^^"'"'^ 

-(1 - x) [3/.ll(xL, 4) + /x- (xL, 4c)] } ■ 
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Figure 8: Size distribution of AlaZr^jSci-^. clusters depending on their number nzr and nsc of 
atoms for a solid solution of concentration = 0.1 at.% and Xg^. = 0.5 at.% at T = 550°C. 
Symbols correspond to KMC simulations and lines to predictions of CNT. 

f^Ah /^Sc /^Zr the chemical potentials in the solid solution: we use a low temperature 
expansion to the second order to calculate them, /x^j, /Xg^. and are the chemical potentials 
in the precipitate. The thermodynamic approach used to calculate them mixes low temper- 
ature expansion so as to correctly describe the majority sublattice which only contains Al 
atoms as well as Bragg- Williams approximation because of the minority sublattice which can 
equally be occupied by Zr or Sc atoms. This approach ensures that for x = or 1, Eq. EH 
leads to the same value for the nucleation free energy as the calculation in the corresponding 
binary alloy fEq. [T3j). 

As the two binary AI^Zt and AI3SC compounds have a close interface free energy and 

(2) 

as the interaction parameter oozv^^ is low in magnitude, we assume a linear interpolation to 
calculate the interface free energy of the AlsZr^Sci.^; compound, 

a{x) = ^^(AlgZr) + (1 - x)ct(A13Sc). (22) 

CNT is not as sensitive as CD to this parameter because not all the thermodynamics is 
based on it. Therefore, one does not need to consider the size dependence of the interface 
free energy as we did for CD. 

The cluster size distribution in the metastable solid solution still obeys Eq. where 
now the cluster free energy is depending on the cluster composition and is given by Eq. 1201 
Comparing this distribution with the one observed during KMC simulations (Fig. |H1), we 
obtain a satisfactory agreement showing that our thermodynamic description of the cluster 
assembly is correct. 

So as to define a classical nucleation rate we use the fact that Sc diffuses faster than Zr. 
We therefore assume that any critical nucleus is growing by absorbing a Sc atom and we 
neglect Zr absorption. A critical size can then be defined for each set of clusters containing 
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Figure 9: Variation with the composition x of the AlaZr^^Sci-a; nucleus of the critical size 
fEq I24|) . the free energy AG* of the critical nucleus, the Zeldovitch factor Z (Eq. OUj) . 
the condensation rate P* (Eq. ITTj) and the nucleation rate J"** (Eq. 1^ for a solid solution of 
concentration x^^. = 0.1 at.% and Xg^ = 0.5 at.% at T = 550°C. 

the same number nzr of Zr atoms. This critical size corresponds to the maximum of the 
cluster free energy (Eq. QUI) with nzr being fixed and is thus given by the condition 

/dAG„Jx)\ ^ 



This leads to a critical size depending on the composition x of the cluster, 

„■(.) . 4.i(.) ^ ( p ::l ;'T'Vr: "^^^f^ t (24) 

We can define now a nucleation rate J**(x) depending on the composition of the nucleus. 
As in the binary alloy, this rate takes the form 

r\x) = N, Z{x) (3*{x) exp {-AG*{x)/kT). (25) 

The Zeldovitch factor Z{x) is defined along the destabilization direction of the nucleus and 
is thus obtained by considering the second derivative of the nucleus free energy with respect 
to TT-sc, keeping nzr constant: 
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Figure 10: Variation of the nucleation rate J'^* (Maximum of J^*{x) given by Eq. l25|) with 
the zirconium concentration for different scandium concentrations of the sohd solution and 
different anneahng temperatures. 

AG*{x) is the cluster free energy (Eq. EUj) of a cluster having the critical size n^{x). As 
for the condensation rate, it is obtained by considering the absorption of Sc atoms by the 
critical nucleus, leading to 

3n^(x) y/^Dsc xl _ ^27) 

271 J C? \ — X 

Writing the nucleation rate in this way, we ensure that when a; = we obtain the same 
expression as that given by CNT in the Al-Sc alloy (Eq. injl . 

The variations of the nucleation rate J**(x) with the composition x of the nucleus are 
shown in Fig. IHl for a chosen supersaturated solid solution. So as to compare results of this 
nucleation model with KMC simulations, we define the nucleation rate J** as the maximum 
of the function J**(x) given by Eq. ESI The critical size used to extract the nucleation 
rate from KMC simulations is the one corresponding to this maximum. When the solid 
solution is highly supersaturated in Sc {x\^ = 0.5 at.% at T = 450 of 550°), our mesoscopic 
models reproduces the variation of the nucleation rate with the zirconium concentration of 
the solid solution (Fig. ^j). The mesoscopic evaluation of J** worsens a little bit when the Zr 
concentration becomes comparable to the Sc one and when the Sc supersaturation becomes 
lower. Nevertheless, the key point is that this extension of CNT manages to predict the 
increase of the nucleation rate with a Zr addition to an Al-Sc alloy, this increase being in 
reasonable agreement with atomic simulations. More advanced models (e.g. the link flux 
analysis and other approximations as described in Ref. [56-58]) are worth to be tried. 

Conclusions 

This study illustrates how a quantitative multiscale modeling of the precipitation kinetics 
can be performed. Using a very limited number of input data, we built for the Al-Zr-Sc 
alloy an atomic model from which mesoscopic quantities like the interface free energy or 



fi\x) 



in 
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the nucleation free energy could be deduced. For the two binary Al-Zr and Al-Sc alloys, it 

was shown that a good agreement can be obtained between the KMC simulations, different 
mesoscopic models (CD and CNT) and experimental data. So as to obtain good predictions 
with CNT, one has to take into account the ordering tendency of the system. In this purpose, 
using low temperature expansions, we propose an analytical expression of the nucleation free 
energy. For the ternary Al-Zr-Sc alloy, we showed that the precipitate inhomogeneity arises 
from the difference between the two solute diffusion coefficients. Although the precipitate 
composition is not known a priori, CNT could be extended to this ternary alloy thanks to 
reasonable assumptions on the kinetics. This mesoscopic model reproduces the increase of 
the nucleation rate associated with a Zr addition to an Al-Sc alloy. 
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